Interaction of vortices in superconductors with k close to \/\/2 
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Using a perturbative approach to the infinitely degenerate Bogomolnyi vortex state for a super- 
conductor with K = l/\/2, T ^ Tc, we calculate the interaction of vortices in a superconductor 
with K close to l/\/2. We find, numerically and analytically, that depending on the material the 
interaction potential between the vortices varies with decreasing k from purely repulsive (as in a 
type-II superconductor) to purely attractive (as in a type-I superconductor) in two different ways: 
either vortices form a bound state and the distance between them changes gradually from infinity to 
zero, or this transition occurs in a discontinuous way as a result of a competition between minima 
at infinity and zero. We study the discontinuous transition between the vortex and Meissner states 
caused by the non-monotonous vortex interaction and calculate the corresponding magnetization 
jump. 

PACS numbers: 74.60.Ec, 74.20.De, 74.55.-|-h, 74.60.-w 



I. INTRODUCTION 

It is widely known that superconducting vortices re- 
pel each other in superconductors of type-II and attract 
each other in superconductors of type-I. The physical ori- 
gin of this phenomenon is the competition between the 
magnetic repulsion of the vortices (dominating in type-II 
superconductors) and the gain in condensation energy of 
overlapping vortex cores producing an attractive inter- 
action (dominating in type-I superconductors). Within 
the Ginzburg-Landau (GL) approximation it can easily 
be seen that the long-range asymptotic behavior of the 
vortex interaction changes its sign at k ~ 1/ \/2. The 
vortex interaction isEJ 



C/int(/) = 2tTc'{k)Ko{1) - ^d{KfKo{V2Kl) 



(1) 



where Kq{1) is the modified Bessel function of zero order, 
and c(k) and d^n) are slowly varying functions of n that 
are equal at k = l/\/2. The detailed profile of ?7int(0 
any / when k goes through 1 / -\/2 was, however, calculated 
only in the 70 's by Jacobs and Rebbia using a special 
symmetry of the z-invariant GL equations at k = 1/ \/2 
discovered by Bogomolnyid. According to Bogomolnyi, 
Jacobs and Rebbi (B JR) the GL energy at k = 1/ \/2 is 
degenerate with respect to any configuration of vortices. 
The sign change of f7int(0 at k = l/\/2 is an exact result 
of the GL theory: at k > 1/^/2 the interaction is purely 
repulsive, at k = 1/V2 vortices do not interact, and at 
K < l/\/2 the iii|t|2raction is purely attractive. 

ExperimcntsQEl, however, reveal a more complex situa- 
tion. The interaction potential L/int(0 close to k = 1/V2 
was sometimes found to be attractive at large distances 
and repulsive at short ones. This non-monotonous pro- 
file of J7int(0 manifests itself as a discontinuous transition 
between vortex and Meissner states and by the existence 



of intermediate mixed-state domains of bound vortices. 
The presence of a local minimum in f/int(0 k ^ l/'s/S 
can be explained (as we do here) by taking into account 
low-temperature corrections to the GL theory, that mod- 
ify the almost flat profile of the intcractionQ. Another 
possibility is to take into account the fluctuations and 
anisotropics in the vortex lattice, that produce an at- 
tractive interaction of the Van der Waals typetl. 

Several ii calculations based on an extendecL, GL 
functionalj were done to clarify this issue. JackobsH cal- 
culated low temperature corrections to the long-range 
vortex interaction (|^) and found that vortices attract 
each other alreadvjin type-II superconductors. Based on 
his results Huberto performed numerical calculations for 
a periodic Abrikosov lattice of vortices and demonstrated 
the non-monotonous behavior of the vortex interaction in 
a vortex lattice. These calculations are consistent with 
numerical variational calculations by BrandttJ based on 
Gorkov equations and solved for vortex lattice configu- 
rarions at all possible values of H, T, and k. Although 
these results reproduce a nontrivial behavior of J7int(0 
at K ~ 1/ v^, it is difficult to survey and interpret them 
in a systematic way because of the cumbersome math- 
ematics and large number of terms in the extended GL 
functional. 

A new method of calculating the properties of supet- 
conductors near k ~ 1/ y/2 has been developed recentlyllil. 
In this approach the degenerate vortex state at k = 1/ V2 
is considered as the starting point of a secular perturba- 
tion theory. A degeneracy-lifting perturbation functional 
in the small parameters 
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is constructed. This approach avoids bulky calculations 
and allows to describe the superconductor with |7|,|i| ^ 



2 



1 in a form that is easy to interpret. 

We will use this perturbation approach to calculate, 
the interaction [/int(0 between two separate vortices and 
between vortices in a lattice, for a superconductor with 
K ~ 1/ when the Ginzburg-Landau theory is extended 
to low temperatures. We study the discontinuous tran- 
sition between the vortex and Meissner states caused by 
the nonmonotonous vortex interaction and calculate the 
corresponding magnetization jump. 



II. PERTURBATION APPROACH 

We start by ou^ning the main elements of our pertur- 
bative approach.U According to BogomolnyiEl, and Ja- 
cobs and RebbiH, in a z-invariant situation, the order pa- 
rameter amplitude [-(/"(r, ri,... ,r„)| of vortices located 
at r = ri, r„ in a superconductor with 7 = 0, t ^ is 
described by the BJR equation 
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V^\n\ilj\=\iPY ^\ + 



(3) 



The magnetic field inside the sample is uniquely related 
to |i/'(r)| via 



6(r) = 1 - IV-Wr 



(4) 



Here, dimensionless variables — /'io, b — y/2KB/Hc 
are used (with '^q the uniform order parameter when the 
external field H = 0). Distances are measured in units 
of the coherence length ^. 

Since at 7 = 0, t ^ the vortices do not interact the 
vortex energy close to this point can be calculated to first 
order in 7 and t by substituting the unperturbed solution 
|V'(r,ri, ...,rn)| of Eq. (||) into the functional 



(5) 



as obtained in Ref. ITT] and defined for the class of Bogo- 
molnyi solutions. Here / = 8tt k'^ T / H'^ + k'^ {T is the GL 
free energy) and hc2 — yf^t^Hcij are the dimensionless 
free energy and the upper critical field. The functional 
accounts for all the terms of the extended GL func- 
tional that are assembled into the terms ji/'l'* and jV'l^ 
with experimentally measurable material coefficients C4 
and cg. The parameter C4 is always positive, whereas cg 
can be both positive and negative. 

To calculate the energy of vortices located at r = 
ri, r„ one should solve first the BJR equation (^) and 
then substitute this solution into the perturbation func- 
tional (|^) . The analytical aspects of this task have been 
discussed in detail in Ref.|ll|; here we solve the BJR equa- 
tion (^) numerically by a finite elements method on an 
adaptive grid by using the Ansatz 



|^(r)|=nk-r,|e 



¥'(r)/2 



(6) 



Using this Ansatz, the BJR equations are reduced to 
a nonlinear Poisson-like equation for Lp and the delta- 
functions are removed. In a weak formulation for the 
finite dimensional space (i^(r) = 'Ylii^i'Pi^y)^ the prob- 
lem is written as 




PiVip X dl - 



(7) 



where, Pj(r) are Lagrange elements on a quadratic gridllj. 
We linearize the exponential function and iterate the re- 
sulting linear part. The dypamic refinement is based on 
the normal Kelly indicatoi£3 for the local error and on 
the error evaluation of the nonlinear exponential part. 

As boundary conditions we use either Neumann (spec- 
ifying or periodic boundary conditions. In the 
Neuman case the normal derivative of (p is taken to be 
d±(p = — 2 X^A: i^-L lii(|r ~ rfel). In the periodic case a 
special constraint is added to the system of linear equa- 
tions for ifi g R, to compensate for the nonperiodicity of 
J2k I'^d'' ~ ^k\)- This algorithm is implemented in C-|— I- 
using the deal-II libraryll3 for the finite elements calcula- 
tions. 



III. TWO VORTEX INTERACTION 

We now use the above method to calculate the in- 
teraction energy L'int(0 between two vortices located at 
ri_2 = ±1/2 by subtracting the self-energy of separated 
vortices 2ei from the two- vortex energy ei^i{l), 



Uint{l) = ei,i(0 - 2£i = eis{l) - £1,1(00). 



(8) 



As follows from (^ the energy J7int(0 '^^'^ be written as 
a superposition of two structure functions Ufc(Z), fc = 4, 6 



C/int(0 = (7 - C4t) U4{1) - cetueil), 



(9) 



which do not depend on the material parameters 7, C4 
and cg. 



Uk{l) 



(1 



|V(r,ri,r2)|'=)-(l-|V(r,ri)|'-) 



-il-\^{r,r2t)]dS. 



(10) 



Here V'(r,ri^2) and ■(/'(r, ri, r2) the one- and two- 
vortex solutions of (|). Note thatllil U2{1) = 0. 

It follows from (^) that the profile of J7int(0 depends 
only on the sign of ce and on the control parameter 



7 ~ 

\c6t\ 



(11) 



In Fig. |i| we show the result of our numerical calcula- 
tions of U4{1) and uq{1). 
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FIG. 1: Structure functions U4,{1) and U(i{l) for two vortices 
separated by a distance I. Appropriate superposition of U4{1) 
and ua(l) gives the vortex interaction energy Uint{l)- The 
distance I is measured in units of the coherence length ^. 

With the numerical method outhned in Sec.^ we can 
calculate the minimum of Uint for 1.4 < d < 2.7. For 
d < 1.4 and for d > 2.7 the minimum cannot be reliably 
found because of the flat profile of C/int(0 small and 
large values of I. To determine the behavior of [/int(0 
I <C 1 and at Z 1 we have used the analytic estimates 
of Ref. |l^ that are summarized below: 

i) the order parameter of the widely separated weakly 
overlapping vortices (1^1) can be approximated as 

|^(r)|^=5?(r + l/2)+.g?(r-l/2)-l, (12) 

where gi (r) is the axially symmetric one-quantum vortex 
solution of the B JR equation. The long-range interaction 
is written as (up to the slow pre-exponential factor u{l)) 

Uintil) = [7 - (C4 + Sce)t] ■ u{l)e-^\ (13) 

and becomes attractive when 7 is smaller than the critical 
value 

7ci = (c4 + dcice)t, del = 3. (14) 

Eq. ( p^ ) has been obtained before in Ref. ^ and general- 
izes Eq. (|l|) to lower temperatures. 

ii) the order parameter of two close-lying vortices with 
almost coinciding cores {I <^1) can be approximated as 

|^(r)| = g2{r) -f \g2{r){\Vf In 52(0, (15) 

where (72(1") is the axially symmetric two-quantum vortex 
solution of the BJR equation. At small I the short-range 
vortex interaction is expanded as 

Ui^t{l) = [0.91(7 - c4) - l.lScet] f + 0{l^). (16) 

It is attractive when the second order term is negative, 
i.e., when 7 is smaller than the critical value 

7c2 = (c4 + dc2Ce)t, dc2 = 1.26. (17) 

We finally compare the energy of a two-quanta vortex 
to the energy of two widely separated one-quantum vor- 
tices. A two-quanta vortex becomes energetically more 
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FIG. 2: The reduced two- vortex interaction energy 
U\i\\.{i) I \c(,t\ at different parameters d = (7 — c/^t') j \c(,t\ (a) 
for C4 > 0, C6 < (crosses indicate the minima) and (b) for 
C4 > 0, C6 > 0. 

favorable when 7 is smaller than the critical value 

7, = (C4 + 4c6)i, 4-T89. (18) 

This value always lies in between 7ci and 7c2. The critical 
parameter dc\ remains unchanged in the case of interact- 
ing multi-quanta vortices whereas the numerical coeffi- 
cients dc2 in Eq. ( p^ ) and dc in Eq. ( p^ ) become larger. 

Combining these numerical and analytical results we 
conclude that depending on the parameters cg and C4 
two scenara are possible for the evolution of ?7int(0 ^ 
function of d: 

I. The case C4 > 0, ce < is shown in Fig. When 
d > del ~ 3 the interaction is purely repulsive. Below dci 
we find an attraction at large distances while the short- 
range interaction remains repulsive. The vortices form 
a bound state with an equilibrium distance lo{d) corre- 
sponding to the minimum of L'int(0- Finally, below the 
second critical value dc2 = 1.26, the short-range inter- 
action becomes attractive too and the vortices combine 
into a two-quanta vortex. 

The distance lo{d) of two bound vortices, shown in Fig. 
^, diverges as — lii{dci — d) for d dd and vanishes as 
y/d — dc2 for d dc2- Special care has to be taken to 
observe these vortex bound states. One can try, e.g., 
to obtain pinned vortex molecules after expelling other 
vortices from the sample by turning off the external field. 
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FIG. 3: Distance between two bound vortices as a function 
of the material parameter d (when C4 > 0, ce < 0). The 
equilibrium vortex distance in the regular triangular vortex 
lattice at discontinuous transition at H — H*i is shown for 
comparison. 



II. The case cg > 0, C4 > is shown in Fig. ||b. The 
order of the critical parameters c?c2, dd is now reversed, 
del < dc2- The vortex interaction now changes its sign 
at negative d. When d > dc2 — —1-26 the interaction 
is purely repulsive. Below dc2 the interaction becomes 
short-range attractive while retaining its long-range re- 
pulsive character. Below dd = —3 the interaction is 
purely attractive. The equilibrium vortex configuration 
in the interval dd < d < dc2 is thus determined by the 
competition between two local minima of J7int(0 at Z = 
and at Z = 00. The formation of a two-quanta vortex be- 
comes more favorable than two isolated and widely sep- 
arated vortices below the critical value dc = —1.89. The 
parameters dd and dc2 then serve as supercooling and 
superheating limits of d. 

The cases considered above exhaust all possible sce- 
naria for the evolution of f7int(0 near the Bogomolnyi 
point. Case I seems to be more realistic given the typ- 
ical values of the material constantalj: C4 ^ 0.1 — 0.5, 
C6 0.2-0. 

We conclude this section by remarking that the vortex- 
antivortex interaction cannot be calculated in this for- 
malism since the vortex- antivortex pair does not be- 
long to the degenerate set of the Bogomolnyi states at 
K = 1/^/2. Moreover, unlike the case of parallel vor- 
tices, both the magnetic and the condensation energy 
contribute equally to the attractive vortex-antivortex in- 
teraction and thus there is no critical point at k = 1/V2. 



IV. VORTEX LATTICES 

The energy of an ensemble of vortices is a nonlinear 
function of the vortex positions. It can be reduced to the 
pairwise vortex interaction (|^) only in the case of large 
vortex separation. Then, in the nearest-neighbor approx- 
imation, the energy of the regular lattice of one-quantum 
vortices with lattice constant ^ ^ 1, calculated with re- 



spect to the uniform Meissner state, can be written as: 

(19) 



b b z 



where ei is the one-vortex energy, z is the coordination 
number, 5/27r is the vortex density. The average magne- 
tization b of the lattice is related to the unit cell area Si 
of the vortex lattice: 



b = 27r/5z. 



(20) 



For equilateral triangular (with z = 6) and square (with 
z = 4) lattices it is: bfj = 27r/P, 6a = iir/P^/S. 

At closer distance between vortices the energy /lat is 
renormalized both because of the next-nearest neighbor 
vortex interaction and because of the nonlinear correc- 
tions. The first contribution can be accounted for by a 
redefinition of C/int(0 as Uint{l) + UintiV^l) for a trian- 
gular and Uint{l) + UintiV^l) for a square lattice. 

However, the nonlinear corrections can only be calcu- 
lated by a vortex-lattice solution of the BJR equation 
(^) . In fact one extend Eq. (|l9|) to the case of arbitrary I 
provided that the contribution (z/2)[/int(0 substituted 
with the interaction energy (2:/2)C/iat(0 calculated from 
the numerical solution of the BJR equation for a peri- 
odic vortex lattice. Similar to the case of two vortices, 
the energy can be expressed as 



{-/-Cit)vi{l)-cetve{l) (21) 
via the lattice-dependent structure functions Vk{l), 



(22) 



v,il)= {l-\Mr)r dSi 



1- |Voc(r)|'=) dS^ 




FIG. 4: Structure functions V4{1) and ve{l) for the square 
(90°), triangular (6O'') and rhombic (45°) lattices with lattice 
constant I. Superposition of W4(0 and veil) gives the lattice 
interaction function {z/2)Uiat{l)- The distance is measured in 
units of the coherence length ^. 

The vortex lattice order parameter |'!/'!(r)| is obtained 
from a numerical solution of Eq. (^), applying the method 
discussed in Sec. ^to a unit cell Si with periodic bound- 
ary conditions. The one- vortex solution ■0 00 (r) | is calcu- 
lated on an infinite unit cell Soc- In Fig. | we show as an 
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FIG. 5: The reduced vortex interaction energy {zUintil))/ 
(2|c6t|) in square (90°), triangular (60"), and rhombic (45*') 
lattices at d = 2.5. 




FIG. 7: Energy / of the different vortex lattices as a function 
of b in an external field ho. The minimum of f(b) corresponds 
to the most stable vortex state at a given ho. The Meissner 
state corresponds to & = 0. The lattice constant / is inversely 

—1/2 

proportional to 6 . The material parameter d is the same 
as in Fig. H. 




FIG. 6: Possible vortex lattice structures considered here, 
f/6 < a/h < I, 7r/f2 < a < tt/2. 



example the structure functions Vk{l) for the triangular, 
square, and rhombic lattices. 

The vortex interaction (z/2)[/iat(0 depends on the con- 
trol parameter d. Here we limit our discussion to the 
most realistic case C4 > 0, cg < 0. We find a non- 
monotonous behavior at d < dd = 3 with long-range 
attraction, short-range repulsion and a minimum at in- 
termediate distances as shown in Fig. || (thus we deal 
with the usual type II superconductor case when d > 3, 
and the type I superconductor when d < 1). We have 
considered all the simple lattices of the form shown in 
Fig. 1^, and, in a non-systematic way, some other lattices 
(multi-quanta, and 4-cluster lattices). Among these, the 
lowest energy has always been found for a triangular lat- 
tice. 

The most important consequence of the non- 
monotonous vortex interaction is the existence of a spe- 
cial class of superconductors that are intermediate be- 
tween type-I and type-II (we call them type I/II). They 
are characterized by a discontinuity of the transition be- 
tween the Meissner and vortex states. A survey of prop- 
erties of these supiesconductors has been given in our pre- 
vious publicationtiJ, where we have demonstrated their 
existence in the interval 1 < d < 3. We discuss now 
the details of the Meissner- vortex transition in type-I/II 
superconductors based on the the expression (|2^) calcu- 
lated above. 



netic field Hq can be investigated by minimizing 

f(b) = /lat(6)"/. = 



b 

2^ 



-Uut{l{b)) - 2TT{ho - h,i) 



(23) 



{fs is the Meissner free energy) over b and different lat- 
tice types. As can be seen in Fig. Q the transition from 
Meissner state with 6 = to the vortex state with finite b 
in type-I/II superconductors occurs discontinuously at a 
critical field h^i that is smaller than the field hd = £i/2tt 
where the penetration an individual vortex becomes en- 
ergetically advantageous. The critical field h^i which 
serves as the lower critical field of the continuous vortex- 
Meissner state transition in type-II superconductors can 
be interpreted as the field of superheating for type-I/II 
superconductors. 

Another important conclusion is that, among the dif- 
ferent types of lattices studied here, the triangular one- 
quantum lattice is the most stable one. We can show 
this numerically for all values of d and ho except for the 
interval 2.7 < d < 3, ho ^ hli, where due to the large 
distance between vortices the accuracy of our calculation 
has been insufficient to draw definitive conclusions. 
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FIG. 8: Jump of magnetization An Am = &(/i*i) at the tran- 
sition between vortex and Meissner state as function of the 
control parameter d (case C4 > 0, cg < 0). 



The stability of the vortex lattice in an external mag- In Fig. ^ we show the magnetization jump 47rAm 
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b{hli) at the h*^ transition as a function of d. Flux 
expulsion varies from almost complete as in type-I su- 
perconductors at d = 1 to vanishingly small at d = 3 
as in type-II superconductors. Surprisingly, we find that 
the jump Am is almost linearly dependent on d between 
these two values. We also show the equilibrium lattice 
constant of the triangular lattice I a = (a/SAz)"^/^ in 
Fig. Fig. 1^ in order to compare it with the distance be- 
tween vortices bound in a pair. At d = 1 the distance 
between vortices in a lattice is minimal and equal to 
= (47r/V3)^/^ ~ 2.69. The vortex distance in the 
lattice is larger than for the vortex pair, because the 
presence of other vortices suppresses superconductivity 
and thus diminishes the pressure of the superconducting 
phase against the vortices. At d ^ dd = 3 the vortex 
distance diverges and the lattice energy can be approxi- 
mated as a superposition of pairwise interaction energies. 
The discontinuity at the h*^ transition can be observed 



experimentally as a spinodal vortex clustering in the ip.- 
termediate/mixed state in thin superconducting platesQ. 
The properties of such clusters are determined by the 
delicate balance between the non-monotonous vortex in- 
teraction C/iat(0 ^'^d the vortex repulsion due to the ad- 
ditional interaction through the stray magnetic field en- 
ergy in the vacuum outside the superconducting plate. 
Detailed calculations of this effect based on the explicit 
form of C/iat(0 currently in a progress. 
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